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Abstract. In this work, we examine effects of large permanent charges 
on ionic flow through ion channels based on a quasi-one-dimensional 
Poisson-Nernst-Planck model. It turns out that large positive perma¬ 
nent charges inhibit the flux of cation as expected, but strikingly, as 
the transmembrane electrochemical potential for anion increases in a 
particular way, the flux of anion decreases. The latter phenomenon 
was observed experimentally but the cause seemed to be unclear. The 
mechanisms for these phenomena are examined with the help of the 
profiles of the ionic concentrations, electric fields and electrochemical 
potentials. The underlying reasons for the near zero flux of cation and 
for the decreasing flux of anion with the increasing of its transmem¬ 
brane electrochemical potential are shown to be significantly different 
over different regions of the permanent charge. Our model is oversimpli¬ 
fied. More structural detail and more correlations between ions can and 
should be included. But the basic finding seems striking and important 
and deserving of further investigation. 


1 Introduction 

Membranes define biological cells. They provide the barrier that separates, and defines 
the inside of a cell from the rest of the world. Membranes are much more than just 
a barrier. They provide pathways for selected molecules to enter and leave cells. The 
barriers must not be perfect or cells would soon die from lack of energy or drown in 
their waste. Biological cells need energy to survive and that is provided (in almost 
all cases) by substances that must cross the membrane. 

Substances cross membranes through proteins specialized for the task. For a very 
long time [20,42] these proteins have been separated into two classes, channels and 
transporters [41], and studied in two traditions, one of electrophysiology [18,47], the 
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other of enzymology [19,37], although the distinction between the two approaches 
was less than clearcut [11]. 

Channels have been viewed fundamentally as charged ‘holes in the wall’ (created 
by the membrane) that could open and close [36] but, once open, the channel followed 
simple laws of electrodiffusion [8]. 

Transporters were viewed as more biological devices, involving conformation 
changes, coupling to energy sources (either ATP hydrolysis or the movement of other 
ions), with quantitative description much more difficult, particularly if the description 
was to be transferable with parameters that were independent of conditions. 

The enormous literature can be sampled in [1-3,37,41]. Structural biology has 
shown that transporters and channels have very similar structures [4,31,38,45]. Bio¬ 
physics has shown that the processes that open and close channels (‘activation’ and 
‘inactivation’) can be coupled to give properties rather like transporters. Physics has 
shown that the electric fields assumed to be constant in classical electrophysiology 
must depend on the distribution of charge in the channel and surrounding solutions 
[9] and so must change with experimental conditions and with location. 

The detailed properties of open channels have not been compared with trans¬ 
porters in the modern literature, as far as we know, certainly not using models 
that satisfy the physical requirement that potential profiles (i.e., electric fields) be 
computed from (and thus be consistent with) all the charges in the system. 

Here we consider a simple model of a permanently open ion channel. (We leave 
gating for later consideration.) Most biologists imagine that if the driving force for 
electrodiffusion is increased-that is to say, if the gradient of electrochemical potential 
across the channel is increased in magnitude - the ffux through the channel should 
increase. We show here that is not always the case. Consider a channel with large 
permanent charge and the ffux of ions with the opposite sign as the permanent charge 
(called counter ions in the ion exchange literature [17] or majority charge carriers in 
the semiconductor literature [34,39,44]. The ffux of ions with the opposite sign as 
the permanent charge in a channel can decrease dramatically as the driving force 
increases - we term this phenomenon as the declining phenomenon. More precisely, if 
the concentration of the ion is held fixed on one side of the channel, and the concen¬ 
tration decreased on the other (‘trans’) side of the channel, the flux of the counter 
ions can decrease if the permanent charge density is large, as we show here. A deple¬ 
tion zone can form that prevents flow even though the driving force increases to large 
values. It is worthwhile to emphasize that if one increases the transmembrane electro¬ 
chemical potential in a different manner, such as, by increasing the transmembrane 
electric potential or the concentration of the ion at one side of the channel, then one 
does not have the declining phenomenon (see Rem. 5.3). 

The decline of flux with trans concentration has been considered a particular, 
even defining properties of transporters, involving conformation changes of state and 
other properties of proteins less well defined (physically) than electrodiffusion [24,43]. 
Declining ffux has been called exchange diffusion (in the early transport literature) 
and self-exchange more recently and is an example of obligatory exchange. Obliga¬ 
tory exchange is a wide spread, nearly universal property of the nearly eight hundred 
transporters known twenty years ago [1,16] with many more known today [40]. Oblig¬ 
atory exchange is ascribed widely to changes in the structure of transport proteins, 
to conformation changes in the spatial distribution of the mass of the protein [14]. 
Obligatory exchange is often thought to be a special property of transporters not 
found in channels. 

The structure of many transporters is now known thanks to the remarkable 
advances of cryo-electron microscopy, recognized in the 2018 Nobel Prize. A trans¬ 
porter (of one amino acid sequence and thus of a perfectly homogeneous molecular 
type) exists in different states. Each state is said to have a different conformation 
meaning, in physical language, that the spatial distribution of mass is different in 


Dynamical Aspects of Mean Field Theories for Electrolytes and Applications 2577 


different states, and the distributions of the different states form disjoint sets, with 
no overlap. The movement of ions is not directly controlled or driven by the con¬ 
formation of mass, however. Rather the distribution of mass produces a distribution 
of steric repulsion forces, and a spatial distribution of electrical forces (because the 
mass is associated with charge, mostly permanent charge of acid and base groups of 
the protein, but also significant polarization charge, as well). It is the conformation 
of these forces that determines the movement of ions. The spatial distribution of 
forces contributes to the potential of mean force reported in simulations of molecular 
dynamics. 

This paper shows that channels with one spatial distribution of mass can have 
properties of self-exchange (for majority charge carrier counter ions) if the density of 
permanent charge is large. The spatial distribution of electrical forces can change and 
create a depletion zone that controls ion movement, while the spatial distribution 
of the mass of the protein does not change. The conformation governing current 
flow is the conformation of the electric field - the depletion zone - more than the 
conformation of mass, in the model considered here. The current flow of counter ions 
is much greater than the current flow of co-ions because there are many more counter 
ions than co-ions near the permanent charge. Transporters almost always allow much 
less current flow than channels. 

It should be emphasized that the depletion zone considered here arises from the 
self-consistent solution of the Poisson-Nernst-Planck (PNP) equations of a specific 
model (large permanent charge, counter ion transport) and that parameters are not 
adjusted in any way to create or modify the phenomena. This is in stark contrast 
to calculations of chemical kinetic models that involve many adjustable parameters, 
without clear physical meaning, and equations that do not conserve current [6 . 

Depletion zones play crucial roles in the behavior of nonlinear semiconductor 
devices although there they usually arise at locations in PN diodes where permanent 
charge (called doping in that literature) changes sign. Depletion zones of the type 
studied here are likely to occur in semiconductors but have received little attention 
because they have less dramatic effects than those in diodes [34,39,44] that follow 
drift diffusion and PNP equations rather like those of open ionic channels [7,9]. The 
possible role of depletion zones in channel function has been the source of speculation 
and experimental verification [10,32]. It is striking that depletion zones can change 
the conformation of the electric field and mimic the obligatory exchange traditionally 
thought to occur only in transporters. Depletion zones can create plastic electric fields 
whose change in shape dominate the flux through a channel of fixed structure. 

Our model is of course oversimplified as are any models, or even simulations 
in apparent atomic detail, of condensed phases. More structural detail and more 
correlations between ions can and should be included. But the basic finding that 
large permanent charge can produce depletion zones and those regions can produce a 
decline of counter ion flux when driving forces increase seems striking and important 
and deserving of further investigation. 

The rest of the paper is organized as follows. In Section 2, we describe the 
quasi-one-dimensional PNP type model and its dimensionless form for ionic flow. 
Assumptions are specified with a key assumption: a dimensionless parameter e defined 
in (7) is small, a dimensionless parameter Q 0 > 0 defined in (8) representing the 
permanent charge is large, and £ is much smaller than 1/Qq. In Section 3, approx¬ 
imation formulas (in small e and v = 1/Qo) for fluxes are provided, which have a 
number of non-trivial consequences. One of the apparent non-trivial consequences is 
that the leading order term Jio of cation flux is zero, independent of the values of 
transmembrane electrochemical potential of the cation (co-ion). The mechanism of 
this distinguished effect of large (positive) permanent charge is examined in details 
in Section 4 with the help of the internal dynamics (profiles of the electric field, 
cation concentrations and electrochemical potential). It turns out the mechanism 
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for Jio = 0 is different over different regions of permanent charge. In Section 5, the 
rather counter-intuitive declining phenomenon - increasing of anion (counter-ion) 
transmembrane electrochemical potential leads to decreasing of anion flux - is shown 
to be consistent with our analytical result. Thus, for the first time (to the best of 
our knowledge), a mechanism of large permanent charge for such a phenomenon is 
revealed. The mechanism is then examined in details again with the help of the 
internal dynamics. We conclude this paper with a general remark in Section 6. 


2 Classical PNP with large (positive) permanent charge 

2.1 A quasi-one-dimensional PNP model 


Our study is based on a quasi-one-dimensional PNP model [12,27,28,33]. For a 
mixture of n ion species, a quasi-one-dimensional PNP model is 


Tl 

®s( £ ' (IWW §) = + «*>)• 

^= 0 , -)A(X)C t ^, k = l, 2 ,...,n ( 1 ) 


where X E [ao, bo] is the coordinate along the axis of the channel and baths, A(X) is 
the cross-sectional area of the channel at the location X, eo is the elementary charge, 
So is the vacuum permittivity, e r (X) is the relative dielectric coefficient, Q(X) is the 
permanent charge density, ks is the Boltzmann constant, T is the absolute tempera¬ 
ture, <f> is the electrical potential, and, for the fcth ion species, Ck is the concentration, 
Zk is the valence, XV (X) is the diffusion coefficient, fik is the electrochemical potential, 
and Jk is the flux density. 

Equipped with system (1), a meaningful boundary condition for ionic flow through 
ion channels (see, [12] for a reasoning) is, for k = 1,2,..., n, 


= V, Ck{ao) = Ck > 0; 3>(6q) = 0, Ck(bo) = IZk > 0. (2) 


Mathematically, we will be interested in solutions of the boundary value problem 
(BVP) (1) and (2). An important measurement for properties of ion channels is the 
I V (current-voltage) relation where the current X depends on the transmembrane 
potential (voltage) V and is given by 


l=Y^ Z sJs(V) 

s —1 


where X/c(V)’s are determined by the BVP (1) and (2) for fixed Ck s and IZk s. Of 
course, the relations of individual fluxes Jk s to V contain more information but it is 
much harder to experimentally measure the individual fluxes JJ s [21]. 


2.1.1 Electroneutrality boundary conditions 

In relation to typical experimental designs, the positions X = ao and X = bo are 
located in the baths separated by the channel. They are the locations of two electrodes 
that are applied to control or drive the ionic flow through the ion channel. Ideally, 
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the experimental designs should not affect the intrinsic ionic flow properties so one 
would like to design the boundary conditions to meet the so-called electroneutrality 

n n 

^z s C s =0 = ^z s 1l s . (3) 

S = 1 S = 1 


The reason for this is that, otherwise, there will be sharp boundary layers which cause 
significant changes (large gradients) of the electric potential and concentrations near 
the boundaries so that a measurement of these values has non-trivial uncertainties. 
One clever design to remedy this potential problem is the “four-electrode-design”: 
two “outer-electrodes” in the baths far away from the ends of the ion channel to 
provide the driving force and two “inner-electrodes” in the baths near the ends of 
the ion channel to measure the electrical potential and the concentrations as the 
“real” boundary conditions for the ionic flow. At the inner electrodes locations, the 
electroneutrality conditions are reasonably satisfied, and hence, the electric poten¬ 
tial and concentrations vary slowly and a measurement of these values would be 
robust. We point out that the geometric singular perturbation framework for PNP 
type models developed in [12,13,22,25-27] can treat the case without the electroneu¬ 
trality assumption; in fact, at the “foot” of the boundary layers the concentrations 
can be determined by the boundary conditions directly and satisfy electroneutrality 
conditions. 


2.1.2 Electrochemical potentials 

The electrochemical potential fik consists of the ideal component /Ajf and the excess 
component fi^ x where the ideal component 


= z fc e 0 $(X) + k B T In ( 4 ) 

Do 

is the point-charge contribution where Co is a characteristic concentration, and the 
excess component fi^ x (X ) accounts for ion size effects. As explained above, although 
not totally physical for ion channel problems in general, we will consider only the 
ideal component in this work that can act as guidance for further studies of more 
accurate models with excess components. 


2.1.3 Permanent charges and channel geometry 

The permanent charge Q(X) is a simplified mathematical model for ion channel 
(protein) structure. It is determined by the spatial distribution of amino acids in the 
channel wall, the acid (negative) and base (positive) side chains, more than anything 
else [9]. We will assume Q(X) is known and take an oversimplified description to 
capture some essence of its effects. For this paper, we take it to be as that in [12], for 
some ao < A < B < bo. 


y\ _ / 0? X E [ao, A) U ( B , bo 
^ X) -\2Qo, Xe(A,B). 



We will be interested in the case where |Qo| is large relative to CA s and TS^’s. 

The cross-sectional area A{X) typically has the property that A{X) is much 
smaller for X E (A, B) (the neck region) than that for X 0 (A, B). It is interesting 
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to note that, in [23], the authors showed that the neck of the channel should be 
“narrow” (small A(X) for X £ (A, B)) and “short” (small B — A) to optimize an 
effect of permanent charge. 


2.1.4 Dielectric coefficient and diffusion coefficients 

We assume that 


e(X) = s r is a constant, and V k (X) = V{X)V k (6) 

for some dimensionless function T>(X) (same for all k) and dimensional constant T> k . 

Note that the assumption V k {X) = V(X)V k is equivalent to the statement that 
T>k(X)/Vj(X) is a constant for k f j. Roughly speaking, the assumption says that, as 
the environment varies from location to location, its influences on the two diffusion 
coefficients V k {X) and V 3 (X) at the same location x are the same; that is, the 
two diffusion coefficients vary from one common environment to another common 
environment in a way so that their ratio is independent of locations. This is not a 
justification of this assumption but only an explanation of what it reflects. 


2.1.5 Main assumptions 

In the sequel, we assume the boundary electroneutrality condition in (3), the form of 
the permanent charge in (5), constant dielectric coefficient and diffusion coefficient 
property in (6), and the electrochemical potential is ideal p k = hfi 1 * n ( 4 )- There is 
a key assumption for this work to be discussed in terms of dimensionless variables 
below. 


2.2 Dimensionless form of the quasi-one-dimensional PNP model 

The following rescaling (see [15]) or its variations have been widely used for 
convenience of mathematical analysis. 

Let Co be a characteristic concentration of the ionic solution. A specific choice 
of Co for the purpose of this paper will be discussed in Remark 2.1 after a key 
assumption. We now make a dimensionless re-scaling of the variables in system ( 1 ) 
as follows. 

2 = £ r £ok B T = X-ap = A(X) = Q(X) 

e l(b 0 - a 0 ) 2 C 0 ’ bo — ao’ (bo — ao) 2 ’ Co 

DO) = V(X), m = ^(.Y), „(,) = _ fjcW ' 

l 

hk(oc) = —-/ifc(I) = z k <fi(x ) + In c k (x). (7) 

k B d 


The dimensionless quantity Q(x) from the permanent charge Q in (5) becomes 


f 0, x € [0, a) U (b, 1] 
\2Q 0) xe(a,b), 


Q(x) 


( 8 ) 
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where 



Qo 


and 


0 < a = 


A 

Vo 


a o 
ao 


<b = 


B 

b~o 


do 

do 


< 1 . 


In terms of the dimensionless quantities, the subinterval (a, b ) C (0,1) corresponds 
to the neck region [A, B]. 

Key assumption. The parameter Qo > 0 is large, the parameter £ is small and is 
much smaller than 1/Qo- 

Remark 2.1. We comment that the assumption is crucial for the mathematical 
treatment in this paper. It is also physically meaningful for the situation we study 
here. Recall that we are interested in the case of large |Qo| A> 1. For definiteness, 
we consider the case Qo > 0- The case where Qo < 0 can be treated in the same 
way. Roughly speaking, if we choose Co in the dimensionless scaling (7) as, say, 

Co = 0(Qq^ 5 ). Then, Qo = 0(Qq^ 5 ) is large, e — O(Q 0 2 ^ 5 ) is small and is much 
smaller than 1/Qo. To have an idea of the size of £ in a physical condition, we 
mentioned that, for example, if bo — do = 25 nm and Co = 10 M, then e nr 3 [5]. 

For fixed Qo? the assumption that £ is small allows one to treat the BVP (9) and 
(10) of the dimensionless problem given below as a singularly perturbed problem. The 
assumption that e is much smaller than 1/Qo insures the validation of an expansion 
of solutions in both parameters £ and 1/Qo- 

We thank one of the referees for pointing out the lack of such a discussion in the 
original submission. □ 

A general geometric framework for analyzing the singularly perturbed BVP of 
PNP type systems has been developed in [12,26,27,30] for classical PNP systems and 
in [22,25,29] for PNP systems with finite ion sizes. 

In terms of the new variables in (7), the BVP (1) and (2) becomes 


e 2 d 


h(pc) dx 
dJ k 


h(x) 


d<p 

dx 


n 


dx 


= 0 , 


S = 1 

J k = D(x)h(x)c k 


^2 z S c s - Q( x ), 

dpk 


dx 


with boundary conditions at x = 0 and x — l 


0(0) = V, Cfc(0) = L k ; 0(1) = 0, Cfc(l) = R k 


where 


T/ e 0 T ^ 

V := -—— V, L k := 


k B T 


Co 


R k := 


n 


k 


Co 


(9) 


( 10 ) 


In this work, we will consider the BVP (9) and (10) for ionic mixtures with one 
cation of valence z\ = 1 and an anion of valence z -2 = —1. It turns out 


a = 


H(a ) J . H(b) 


HO) 


and [3 = 


HO) 


where H(pc) = / 

Jo 


X 


1 


o D(s)h(s) 


ds 


are the key parameters that characterize the effect of channel geometry and location 
of permanent charge. The physical meanings of these parameters could be seen from 
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the special case when h(x) = ho and D(x) = Do are constants. In this case, 


H{x) = 


x 


Doho 


which is proportional to the (scaled) length x of the region over [0, x\ of the channel 
(conductive material), inversely proportional to the (scaled) cross-sectional area ho of 
the channel and to the (scaled) diffusion coefficient or electric conductivity Do] that 
is, H(x) represents the resistance of the portion of the channel over [0,x]. 


3 Approximations of fluxes [46] 

Before presenting the approximations of fluxes obtained in [46], we briefly describe 
the geometric or dynamical system approach for BVP of (9) and (10) established in 
[12] for n = 2 (see, [27] for general n where a treatment of both the limiting fast 
(inner) system and the limiting slow (outer) system is fully developed); in particular, 
we briefly comment on nature of Debye layers of the solution and refer the readers 
to the above-mentioned references for details. 

Introduce u = £<f> and w = x, system (9) can be recast as, for k = 1, 2, 


1 ■ r>r ^ £h'(w) 

£(f> = U, £U = —Z\C\ — Z2C2 — Q(W) -——— U, 


£C-k — 


&Jk 


D(w)h(w) 


h(w) 
Jk = 0, w = 1, 



where the symbol dot denotes the derivative with respect to the x-variable. The 
equation w = 1 is augmented so that system (11) is autonomous and can be treated 
as a dynamical system with phase space M 7 and state variables (</>, u , ci, C 2 , Ji, J 2 , w). 
The BVP is then reduced to a connecting problem: finding an orbit of (11) from 


Bo = {(V,^,Li,L 2 ,Ji,J 2 ,0) : arbitrary u, Ji, J 2 } 


to 


B\ — {(0, u , i?i, i? 2 , Ji, J 2 ,1) : arbitrary u, Ji, J 2 }. 


An orbit for the connecting problem will be called a connecting orbit of the BVP. 
Note that, in constructing a connecting orbit of the BVP, one needs not to track the 
x -variable since w = x is encoded in the orbit of system (11). In particular, whenever 
a connecting orbit is constructed, the variable w = x goes automatically from 0 to 1. 
A great advantage of the connecting problem is that one can scale the vector field by 
arbitrary positive functions which could depend on the state variables and the phase 
portrait stays the same. This is one of the two critical ingredients in the geometric 
framework established in [27] specifically for PNP. 

System (11) is called the slow (outer) system of the singular perturbed system. 
The fast (inner) system is, for k = 1,2, 


x , / \ £h’(w) 

(p = U, U = ~Z\C\ — Z2C2 — Q\UJ) -— —U. 


C 


k 






h(w) 


j’k = 0, w' = £, 


D(w)h{w) ’ 


( 12 ) 
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Fig. 1. An illustration of a singular connecting orbit projected to the space of (u,z\C\ + 
Z 2 C 2 , w). Boundary layers at w = 0 and w = 1 exist if electroneutrality boundary conditions 
are not assumed. 


where prime is the derivative with respect to the fast variable £ = xje. For e > 0, the 
fast system (12) has the same phase portrait as that of slow system (11). 

A singular (zeroth order in e) orbit of the connecting problem will be a union of 
slow orbits (outer solutions) for the e = 0 limiting slow system and fast orbits (inner 
solutions) for £ = 0 limiting fast system. 

In view of the jumps of permanent charge Q(x) in (8) at x = a and x = b, the 
construction of singular orbits is split into three intervals [0, a], [a, b], [6,1] as follows. 
One preassigns six (unknown) values of (</>, Ci,c 2 ) at x = a and x = b: 


H a ) = 0 a , ci(a) = c“, c 2 (a) = eg; <j>(b) = (j) b , ci(b) = c\, c 2 (a) = eg 


(13) 


Then these values determine two (boundary) sets at x = a and x = b as 


B a = {(0 a ,u,Ci,C2 , Ji, J 2 ,a) • arbitrary u,J ll J 2 } 1 


and 


Bb = {(0 D , c?,4, Ji, J 2 , b) : arbitrary u, Ji, J 2 }. 


On each subinterval, a singular orbit can be constructed (see [12]) and typically 
consists of two singular (inner) layers and one regular (outer) layer (see Fig. 1 for an 
illustration). 


(i) On interval [0, a], a singular orbit from Bq to B a consists of two singular layers 
located at x = 0 and x = a, denoted as Tg and T l a , and one regular layer A/. 
Furthermore, with the preassigned values </> a , c\ and the flux J\. and ui(a) 
are uniquely determined so that 


{(j) a ,ui{a),c a 1 ,c a 2 , j[,J l 2 ,a) e B a . 
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(ii) On interval [a, b], a singular orbit from B a to B 5 consists of two singular layers 
located at x = a and x = b, denoted as T r a and T l bl and one regular layer A m . 
Furthermore, with the preassigned values (</> a , cj, C 2 ) and (</> 5 , c\, c^, b), the flux 
J™, u m (a) and u m (b) are uniquely determined so that 


G £? a and e B b . 


(iii) On interval [b, 1], a singular orbit from £>5 to T>i consists of two singular layers 
are located at x = b and x = 1 , denoted as T£ and T^, and one regular layer A r . 
Furthermore, with the preassigned values </> 5 , c\ and <4, the flux J k and ?z r (b) 
are uniquely determined so that 


The matching conditions for singular orbits of the full connecting problem are 


J l k = J k l = J k for k = 1 , 2 , iq(a) = u m (a) and u m (b) = u r (b). 


(14) 


The number of conditions in (14) is six, which is exactly the same number of unknowns 
preassigned in (13). In this way the singular connecting problem is reduced to the 
governing system (14) (see [12] for an explicit form of the governing system). 

Once a singular orbit is constructed and with some checkable transversality 
conditions, one can justify the validation of the singular orbit for £ > 0 small [12 
It should be pointed out that, typically, the internal layer T^ U T^ at x = a 


tias 


a Conner at the joint so it is a double-layer (see Fig. 1). Similarly, the internal layer 
U T^ at x = b is also a double-layer. 

We now recall the results on approximations of (</>, Ci, c 2 , Ji, J 2 ) and (i k s from 
46] for the case where z\ = 1 and Z 2 = — 1 with L\ = L ) = L and B\ = /? 2 = R- 
We comment that, in [46], the electroneutrality boundary conditions L\ = L 2 and 
R\ — R 2 are not required. 

Based on the assumption that Qo > 0 is large, £ > 0 is small and is smaller than 
1/Qo? one can expand fluxes in £ and in v = l/Qo near £ = 0 and v = 0. 

First, one expands the fluxes in £ as 


Ji(£, v) = J\[y) + 0(e) and J 2 (e, v) = hiy) + O(e), 

where Jfc(^), depending also on the parameters (V, L, R, H(l), a, /?), are the zeroth 
order in e terms of the fluxes. Then, one expands Jk(v) about v as 


J\(y) = J 10 + J 11 V + 0(iz 2 ) and J 2 (z/) = J 20 + J 21 V + 0(v 2 ). 

Thus, 

Ji{e, v) = J 10 + J\\v + 0[y 2 ,e) and J 2 (e, v) = J 20 + J<i\v + 0(y 2 , e). (15) 

Note that J 10 and J 20 , depending on (V, L, i?, if(l), a,/3), contain the leading order 
effect of the small v (or equivalently, large Q 0 ). 
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The following result about the fluxes is established in [46] 

Proposition 3.1. One has 


J\o — 0 , 


Jn = 


1 


J 20 — 


J21 — 


2iT(l)08 
2 \J LR 


(1 — (3)L + aR 
ol) \(1 - (3)V^L + aVR 

1 


(e v L-Rf 


H( 1 ) (1 - f3)VL + aVe~ v R 
((3 - a)e v LR((l - f)L + aR) 
H(l) ((1 — f3)Ve^L + ay^R) 3 


(Ve~ v L-VR), 


(Ve~ v L-VR) 


+ 


(e v L — R)( — V + InL — In R) ((1 ~/3)L + afl) 3 
4(/3 - ct)ff(l)(\4rU - W?)((l - P)\feVL + ay/R) 

e v L — R ( (1 -P)L + aR 

2(8 - a)H(l) In- BWe v L + 




(16) 


A distinct implication of J10 = 0 in Proposition 3.1 is that large (positive) per¬ 
manent charge Qq (or small v = 1/Qo) inhibits the cation flux. We will provide a 
detailed discussion in Section 4 on what happens to the internal dynamics that is 
consistent with this conclusion. 

To this end, we recall another immediate consequence of (16) (see [46] for more). 

Corollary 3.2. [Current Saturation] For large permanent charge Qq (small v = 
1/Qo) an d to the leading order terms in v, individual fluxes Jk ; s, and hence, the total 
current I saturate in \V\; more precisely, one has J20 is decreasing in V, concave 
downward for V < Vf and concave upward for V > Vf for some Vf, and Jn is 
increasing in V, concave upward for V < V* and concave downward for V > V* for 
some V*. Furthermore, | J201 an d | Jn| are uniformly bounded in V with a bound that 
can be determined from the limits 


lim J 20 = 

V —^+00 


lim Jn = 

V —^+00 


1 


2 R 


i- pH(iy 


lim J 2 i = 

V— ^+cx) 


r 12 L 

lim J20 = -T77TT, 

V—y—00 a H( 1) 

1 ((1 -P)L + aR)' 


lim Jn = — lim J? ^ = 


V —y — 00 


V —y — 00 


21 


(I-PY 2H(1)((3 — a) 
1 ((1 -f3)L + aR) 2 
2H(l)([3-a) ' 


Proof. All the above statements can be derived from (16) easily. 

Note that it is NOT the case if the permanent charge is small [23] 


(17) 


□ 


Remark 3.3. We want to emphasize that the formulas in (16) are established for 
bounded V. Strictly speaking, the formulas are not justified as \V\ —>■ 00. The limits 
are presented, together with the monotone properties of J20 and Jn in V, to allow the 
conclusion that, for bounded V, J20 and Jn are bounded with a bound independent 
ofV. □ 


We end this part by comparing our analytic results with the numerical results 
in pages 180-185 of [35] (to be simply referred to as [35] below). In [35], the author 
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presented a numerical study of the BVP with more or less the same setup as in our 
work: Dirichlet boundary conditions and the same form of the permanent charges, 
except that the nonzero value of the permanent charge, denoted there by — N (so 
that —N = 2Qo), was taken to be negative (N > 0). In particular, the co-ion in [35] 
is the negatively charged ions and, in our case where Q o > 0, the co-ion is the cation 
(label with subscript 1). The boundary concentrations are special and fixed in [35] 
and are taken to be L = R = 1 (in terms of our notation). Numerical results were 
presented for relevant quantities as £ varies, as V varies, and as N varies. 

It turns out our analytical results are consistent with the numerical results, say, 
reported in Figure 5.3.2 b, c of [35]. (We only examined the zeroth order terms in £ 
so cannot compare with the numerical I-V curves with different £ values reported 
in Fig. 5.3.2a.) More precisely, our analytical results for the approximated current I 
would be 


I = J\iv ~ (J 20 + J 21 V) up to 0(z/ 2 ,e), 

and it can be checked from our formulas for J20, J 21 and J 11 with L = R = 1 fixed 
in the numerics in [35], that \I\ is decreasing in N = —\/v for fixed V, which is 
consistent with the qualitative behavior reported in Figure 5.3.2b. 

Figure 5.3.2c in [35] reports numerical results for the co-ion (cation in our case) 
transport number t\ = |Ji|/|/| with the values N = 10, 20, 40. These values are 
moderate but could be viewed as large. In our case, for large Qq = 1/v, the cation 
flux and anion flux are approximated as 


Ji — J10 T Juv H~ 0(^ 2 ,<s) — J\\v + 0(z^ 2 ,<s), J 2 — J20 T J 21 V ~b 0{v 2 ,e). 


Thus, 


Tl 


Jl 

I/I 


Juts T d)(z^ 2 , s) 


— J 20 — J 21 V 4~ J\\v T 0(^ 2 ,s) 


JiM 

IJ20I 


+ 0(is 2 , 


Note also that, for L = R = 1 and h(x) = 1, one has, from (16), 


J 11 = 


k 


( e v /2 + 1)2 


(e v - 1) . j20 _ - 


e V/2 _|_ l 


1 + e W 2 


e )- 

1) 


where kj 1 s are constants results from a, /3 and H( 1) in our formulas and are 
independent of V. Thus, for V > 0 but bounded, 

\Jnv\ k\ (e v / 2 — 1)(1 + e~ v / 2 ) 2 2k\ 

Tl “ P20I = F (1 — e~ v / 2 ){e v / 2 + 1)~N = k^N’ 

which is consistent with the numerical results presented in Figure 5.3.2c. The number 
Ti is nearly constant in V and, for fixed V, it decreases as N increases. Also the trend 
in the figure indicates that r\ is close to zero as N becomes large. 


4 Internal dynamics for J 10 — 0 


It follows from the Nernst-Planck equation in (9) that 



1 


Ji 


0 


D \(x)h(x)c\(x) 


dx = fi i(0) — (i i(l). 















Dynamical Aspects of Mean Field Theories for Electrolytes and Applications 2587 


Thus, J\ has the same sign as that of fl i(0) —/2i(l); in particular, if /xi (0) — /ii (1) ^ 0, 
then Ji 0. 

For the setup of this paper, there are three regions of permanent charge Q(x): 
Q(x) =0 for x G [0, a) U (6,1] and Q(x) = 2Qo for x G [a, b] with large Q o or small 
v = 1/Qo- A major consequence of large Qo is, from Proposition 3.1, up to the leading 
order in v = 1/Qo and in £, the flux of cation is Jio = 0, even if the transmembrane 
potential ft i(0) — fl i(l) ^ 0. We will reveal the internal dynamics that lead to this 
conclusion. To do so, we will discuss what happens over each subinterval based on 
the approximated (zeroth order in e) profiles. Let 

(<p(x; e, v), c k (x; e, v), J k (e, v)) = (<p(x; v), c k (x; v), J k (v)) + O(e) (18) 

be the solution of the BVP (9) and (10). For v > 0 small, one has the following 
expansions 


Ji{v) — Jio + A 0(v 2 ), 

where Jio = 0 and Jn are given in (16). We will also provide figures for the profiles of 
rescaled electrical potential <f(x; v), cation concentration ci(x; v) and electrochemical 
potential pi(x]u) of the cation, respectively. The parameter values used are 

e 0 =1.6022 x 10- 19 (C), k B = 1.3806 x 10“ 23 (JR -1 ), T = 273.16 (K), 

V =0.01 (V), C = 10 (M), U = 2 x 10“ 4 (M), Qo = 2 x 10 3 (M), C 0 = 20 (M), 
ao =0, b{) = 25 (nm), A = (6 0 — ao)/3 (nm), B = (6 0 — a o)/2 (nm). 

We comment that the value for Qo is not physical. We use this value for some con¬ 
sistence of the size requirement of £ and Qo in the key assumption given after the 
rescaling (7), which is required for a rigorous expansion used in this work. The com¬ 
parison of our analytic formulas to some numerical results in [35] at the end of the 
previous section indicates the analytical results could be valid for a more realistic 
range of the parameters. The corresponding dimensionless parameter values are 

V = 0.425, L = 0.5, R = 10" 5 , a =-,b=-, v = 10" 2 , £ « 2 x 10“ 3 (19) 

3 2 

The scaled cross-sectional area h(x) is taken to be 

{ 7r(— x T ro T a) 2 , x G (0, a) 

7rrQ, x G (a, b) 

tt(x + ro — b ) 2 , xG (&, 1). 

The same choice of the parameter values will be used for figures in Section 5. 

Notice that the whole interval is (0,1), and it is divided into three subintervals 
(0, a), (a, b) and (6,1). Permanent charge is located over (a, b). The radius of the neck 
of ion channel is ro- In the figures, we let ro = 0.5. 

4.1 Internal dynamics over the interval (0, a) 

The leading order terms of (0,Ci,/2i) in (18) are derived in [46]. One has 

Proposition 4.1. For x G (0,a), 


(j){x] v) = (j)o(x) A 0(u), Ci(x] v) = Cio(x) + 0(y), [i\[x\ v) = /i 10 (x) A O(iz), 
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where 

c 10 (x) = L - H(x), <po(x) = V - In Cl0 ( x \ 

P-io(x) = <po(x) + In c w (x) = V + In L. 

The following direct consequence explains the reason for Jiq = 0 over (0,a). 

Corollary 4.2. Over the interval (0,a), cio(x) = 0(1) but /2' 10 (x) = 0, and hence, 
one has Jiq = 0. 


4.2 Internal dynamics over the interval (a, b) 

It follows again from [46] that one has the following approximations. 

Proposition 4.3. For x E (a, b), 

(j)(x; v) = —In v + 4>o(x) + 0(v), 

Ci{x\v) =cio(a;) + cu(x)u + 0(u 2 ), 
pi(x;v) =pio(x) + Op), 


where 


, , , , 2 e v L . , n 

<Po{x) = In i cio(x) = 0, 


L o 


hio(x) = <f> o(x) + lncn(x) = In 


Clip) = \aI-J 11 {H{x)-H{o)), 

H(b)-H(x) v H(x)-H(a) \ 

H(b) - H(a) H(b) - H(a) ) 


with 


\/eVL((l ~^)L + afl) 

° ~~ (1 - j3)VerL + as/R ’ 

Note that the expansion of <f(x] v) in v is not regular with the term — Inn over 
this interval and the zeroth order term pio(x) involves the first order term Cu(x). 

One has the following immediate consequence providing a mechanism for Jiq = 0 
over (a, b) that is different from that over (0, a). 

Corollary 4.4. Over the interval (a, b), c\q(x) = 0, and hence, Jiq = 0. Furthermore, 
pi io(a) = V + InL and pio(b) = In R, and, as R decreases, piQ(a) — flio(b) = V + 
InL — Ini? increases (without an upper bound independent of R). 

Proof. It follows from Proposition 4.3 directly. □ 


4.3 Internal dynamics over the interval (5, 1) 
Proposition 4.5. For x E (b , 1), 


<t>(x] v) = <t>o(x) + O(iz), c i(x; v) = Cio(x) + 0(y), [i\[x\ v) = pbi 0 (x) + 0(n), 
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where 


cioO) =R + — - H(x)), (j) 0 (x) = In R — lnci 0 (x), 

Jio (x) = (j) o (x) + In cio {x) = In R. 

The following consequence shows a mechanism for Jio = 0 over ( 6 ,1) that is 
similar to that over ( 0 ,a). 

Corollary 4.6. Over the interval (6,1), cio(x) = 0(1) and j' 10 (£) = 0, and hence, 
Jio = 0 . 


4.4 Summary of mechanism for J 10 = 0 

From the above discussion, we conclude that the mechanisms for Jio = 0 are different 
over each subintervals. More precisely, one has, 

- over the first interval (0, a), Jio = 0 is the result of constant ftio(x) = ft io(0) 
(so that fiio(x) = 0) while Cio(x) f 0; 

- over the last interval (6,1), Jio = 0 is the result of constant ftio(x) = ft io(l) (so 
that ftio(x) = 0) while c\q{x) f 0; 

- over the interval (a, b) where permanent charge is large, Jio = 0 is the result of 
cio(x) = 0 while ft\o(x) is not a constant, in particular, the drop of ftio(x) over 
(a, b) equals its drop over (0,1), in fact, ftio(a) = fti(0) and Jio (b) = Ji(l). 

Here we provide the profiles of concentration ci(x; v) (Fig. 2), electrical potential 
cj)(x; v) (Fig. 3) and electrochemical potential fti(x; v) (Fig. 4) over the whole interval 
[0,1]. The concentration Ci(x; v) and the electrical potential v) are not continuous 
at x = a = l/3 and £ = 5 = 1/2 due to the discontinuity of the permanent charge at 
these point. On the other hand, the electrochemical potential fti(x; v) is continuous. 


5 Declining phenomenon and internal dynamics 

In this section, we will show that large permanent charge is a mechanism for the 
declining phenomenon described in the introduction. Recall that, by the declining 
phenomenon , we mean the following: 

For fixed V and L, as R decreases, the flux of counterion (J 2 in the setting since 
Qo > 0 / decreases monotonically without a positive lower bound. 

Remark 5.1. The phenomenon was well-known in the physiology community. Unfor¬ 
tunately, we could not find references stating precisely this phenomenon. We have 
contacted many leading experts who are all recognizing this phenomenon. Some 
experts mention this phenomena as an example of ‘exchange diffusion’ and/or long 
channel phenomena. □ 

This phenomenon is rather counterintuitive since the (scaled) transmembrane elec¬ 
trochemical potential for the counter-ion p 2 ( 0 ) — ^ 2 ( 1 ) = —V + InL — In R increases 
(without bound) as R decreases. 
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Concentration of Cation in (0,1) 

0.6 - 1 - 1 - 1 - 1 - 

0.5 -- 

0.4 - 

0.3 - 

2 

o 

0.2 - 

0.1 - 

0 - - - 

- 0.1 - 1 - 1 - 1 - 1 - 

0 0.2 0.4 0.6 0.8 1 

x 

Fig. 2 . Profile of ci(x; v) over [0,1]. Note that c\q(x) — 0 for x £ (a, b) and, for R — 10 5 
small chosen for the numerics and since J 20 = 0(y/R), cio(x) ~ L = 0.5 for x £ (0,a) and 
cio(x) = o(Vr) is close to zero for x £ (b, 1). 


Electropotential in (0,1) 



Fig. 3. Profile of <j>(x] v) over [0,1]: For v — 10 2 and R = 10 5 small chosen for the 
numerics, <j)(x\v) ~ V = 0.425 for x £ (0,a), (f)(x;v) = 0(Inv) for x £ ( a,b ), and <j)(x\v) = 
0(lni?) for x £ (6,1). 
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Electrochemical Potential of Cation in (0,1) 



Fig. 4. Profile of fii(x;is) over [0,1]. For v — 10 2 and R = 10 5 small chosen for the 
numerics, ~ F + In L — —0.214 for x E (0, a) and JX\(x\v) ~ In R= —11.5 for x E 

(5,1). Note the large drop of jxi over the interval (a, b) = (1/3,1/2) and the nearly constant 
values over the other two subintervals. 


5.1 Experimental phenomena are consistent with our analysis 


The formula for J 2 o in (16) actually justifies the declining phenomenon, up to the 
leading order J 2 o in v for small v (or for large Q o); that is, 

Proposition 5.2. As a function of R, the leading order term J 2 o of the flux is 
monotonically decreasing and concave downward. Furthermore, as R deceases, |J 2 o| 
deceases (without lower bound independent of R). 

Proof Indeed, from the expression of J 2 o in (16) and treating J 2 o as a function of w 
where R = w 2 for convenience, one has 


^ 20 <» = 


\[L V e v Lw — w 2 


1 ) (1 — fl)y/L + aVe v 


w 


It is clear that J 2 o(^) 0 + as w —0 + . Note that the derivative of J 2 o in w is 


^2o( w ) ~ 


1 (1 - /3)LVe^ - 2(1 - fl)y/L‘ 


w — aVe v w 2 


2H(1) 


[(1 - (3)\/1j + aVe v 


w 


Thus, from the expression of the numerator, if w is smaller than some wq, then 
^2C)( W ) > 0, and hence, as w —> 0 + (or equivalently, R —> 0) monotonically over the 
interval [ 0 ,rco], J 2 o(^) —> 0 monotonically. 

It is not hard to show that for R > 0 but smaller than some positive value, the 
graph of J 20 as a function of R is concave downward. □ 
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In Figure 5, the horizontal axis is for R and the vertical for J 20 . We fix V = 
0.425 and L = 0.5 as in (19), and vary R £ (0,10~ 3 ]. The monotonicity and concave 
downward features of the graph are apparent. 


Remark 5.3. We comment that, for large permanent charge, the declining curve 
phenomenon occurs when the transmembrane electrochemical potential /^(O) — ^ 2 ( 1 ) 
is increasing in a particular way ; that is, as R deceases with fixed V and L. If one 
increases the transmembrane electrochemical potential /^(O) — fl 2 ( 1 ) in a different 
manner, for example, as | W| increases or as L increases, Corollary 3.2 shows that the 
declining curve phenomenon does not happen. 

It is also important to note that, when the next order term J 21 V is considered, 
then, as R decreases, |J 21 1 increases. But, if R and z/ln R are small, then \J 21 v 
stays small. Thus, only when v is very small (Q 0 is very large), is the term J 21 V not 
significant, and hence, the term J 20 dominates the described behavior. □ 


5.2 Mechanism of declining phenomena from the profiles 

Recall, from the Nernst-Planck equation in (9) that 

d 

— J 2 = D 2 (x)h(x)c 2 (x ; v)—ii 2 {x] v 
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Since D 2 (x) and h(x) are fixed, we will treat them as of order 0(1) quantities so that 
they do not contribute much to the near zero flux scenario when R is small. Thus, as 
far as the near zero flux mechanism is concerned, one has 


-J 2 c 2 (x ; v)ii f 2 (x\ v) 


( 20 ) 


One sees that the gradient \F 2 [x\v) of the electrochemical potential is the main 
driving force for the flux J 2 . Intuitively, large drop of (or transmembrane) electro¬ 
chemical potential ft 2 ( 0 ) — p 2 (l) of jH 2 produces large flux J 2 . In this sense, the 
declining curve phenomenon is rather counterintuitive. A careful look at (20) reveals 
that there is only one possibility for the declining curve phenomenon; that is, when¬ 
ever p! 2 {x\ v) is large, c 2 (x ; v) has to be much smaller in order to produce a small flux 

J 2 \. We will apply the analytical results of the internal dynamics from [46] to show 
that this is indeed the case. 

For fixed V and L, jH 2 (0) — ft 2 ( 1 ) = — V + InL — Ini? « — In R 1 for small R. 
We need to understand 

(i) HOW the electrochemical potential jH 2 drops an order O(-lnR) 1 over the 
interval x E [ 0 , 1 ]; 

(ii) HOW J 2 can be small, for small v and R , at every x E [0.1]; 

(iii) Most importantly, HOW the above two things, with the constraint (20), can 
happen simultaneously. 

There are two small parameters v and R in this consideration. The relative sizes 
of these two parameters is relevant for J 2 q being a good approximation for J 2 . It 
turns out we will need vhiR ^ 1. In order to show this, we will examine the profiles 
of (cf)(x] v), c 2 (x ; z/), fl 2 (x- v)) up to order 0( v). 

We now discuss what happens over each subinterval based on the approximated 
(of zeroth order in e) functions of profiles. To do so, let 

£, v),c 2 (x\ e, v), J 2 {e, 0) = iy )> c 2 (x; v), V(^)) + 0(e) 


be the solution of the boundary value problem. For v > 0 small, one has the following 
expansions 


Jziv) — j 2 0 + ^21^ + 0(z^ 2 ), 

where J 20 and J 21 are given in (16). The expansions in v for <\>{x\ v) and c 2 (x; v) are 
not regular oyer (a, b ) and are qualitatively different over the subintervals (0, a), (a, b ) 
and (5,1). They will be given explicitly in each subsection below for us to understand 
what happens over each subinterval. 


5.2.1 Internal dynamics over the interval (0, a) 

The leading order terms of (cj),c 2 ) are derived in [46]. One has 

Proposition 5.4. For x E (0,a), 

4>(x; v) = 4>o(x) + <j>i{x)v + 0(u 2 ), c 2 (x; v) = c 20 {x) + c 21 (x)v + 0(u 2 ), 
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Concentration of Anion in (0,a) Electropotential in (0,a) 



Electrochemical Potential of Anion in (0,a) 



Fig. 6. Profiles of c-i(x\ //), 4>(x\ //). and /.(■>(.:/;; v) over interval [0, a]. For v = 10 2 and 
R = 10 -5 small used for numerics, ci(x\v) = L + 0(y/R) ~ 0.5, <f>(x;n) = V + 0(\/R) ~ 
0.43, and p, 2 (x ; v) ~ — V + In L = —1.12. 


where 


Mx) =v 


C 20 (x) =L 



<t> i(x) 

C2l(x) 


C 2 i(x) 2 Jn C 20 (a:) 

c 2 o(x) J 20 n i 5 


•A i + J21 
2 




In particular, ^(x; v) = ^20 0*0 + P 2 i(x)v + 0{u 2 ), where 


H20(x) 


— V + 2 In C 2 o(;c) — Ini, 


A21P) = 2 


/ c 2 i(x) 

V c 2 o(®) 


Jn 

J20 



Note that C 2 o(x) = cio(x) and C 2 i(x) = cn(x) for x E (0,a). 

Figure 6 shows profiles of 02 ( 2 ;; v), <f>(x; v) and p2(x; v) over the interval (0, a). 

Corollary 5.5. Over the interval (0, a), C 2 (x;v) = 0(1) but 


/40;z/) 


J20 

h(x)c 2 o(x) 


+ 0{y In R) = 0(VR, v In R). 


Therefore, from (20), 

d 

J2 ~ C2(x; v)—jl2{x] v) = 0 ( 1 ) 0 (VR , v In R ) = 0 (VR , z/ In i?) 

dx 
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and fi 2 (x; v) drops an order of 0(\fR, v In R) over the interval (0, a). 

In particular, for zeroth order in n, J 2 q = 0(VR) and p 2 q(x) = 0(\fR) over the 
interval (0, a). 


5.2.2 Internal dynamics over the interval (a, b) 


It follows from [46] that 

Proposition 5.6. For x E (a,b), 

(j>{x\ v) = -Inn + </>o(#) + (fi(x)u + 0(v 2 ), 
c 2 (x;is) = ^ + - Jn(H(x) - H{a)))v + 0{v 2 ), 

where Aq is given in Proposition 4-3, and 


<p 0 {x) 


1 2eVL A ( \ 

In— <j>i(x) 


a 


A 0 + h^W(x)-H(a)), 


l o 


A i = 


a(j3 — a)e LR((1 — j3)L + aR) , ^Z\FJ _ 

2((1 - /3)VVU + ay/R) 3 
a(V - In L + In R){{ 1 -/?)£, + aR) 3 (e v L - R) 


+ 


8{/3 - a){'/e =v L - VR )((1 - /3)VWL + aVR) 


3 ’ 


0 


a 


In 


At 


In 


e v LR (P- a ( L ^p) W ~ P)A ZL + gy/R A x 


L 


a 


+ 2 


e v R 

(1-(3)V^ V LA 1 (3 


a 


Vr 


A 


0 


a 


y/~R Aq ex 


(L — Aq) + L. 


In particular, p 2 (x; v) = p 2 o(x) + p 2 i(x)u + 0(n 2 ) where 


A 2 

h 20 (x) = -0o (x) + In 2 = In , p 21 (x) 


~ 01 + Aq — 



(H(x)-H(a)). 


The profiles of c 2 (x; n), <j>(x; n), and p 2 (x ; v) are shown in Figure 7. 

Note that, over this interval, both <f(x; v) and c 2 (x ; v) are singular in za To under¬ 
stand what happens to the internal dynamics for anions over the interval (a, 5), we 
need to resolve this singularity by considering at least 0(n) -terms. 

Corollary 5.7. Over the interval (a,b), c 2 (x;n) = 0(1/v) 1 and p 2 (x;n) = 

O(lnR) 1 but 


r 2 {x-,v) ss /j,' 21 (x)u = - 


J 20 

2 h(x) 


0{vXR). 


Therefore, from (20), 

J 2 ze 0(l/v)0(vyf~R) = 0 (VR) and p 2 (x;n) only drops 0(v\fR) over ( a,b ). 

Remark 5.8. Note that the drop of p 2 (x) over the interval (a, b) is much less than its 
drop over the interval (0, a). But both drops are small and contribute nearly nothing 
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Concentration of Anion in (a,b) Electropotential in (a,b) 



Electrochemical Potential of Anion in (a,b) 



Fig. 7. Profiles of C 2 (x;z/), </>(x;z/), and /i 2 (x over interval [a, b]. For v — 10 2 and 
R = 10 -5 small used for numerics, C 2 (x; v) — 2/v + 0(y) ~ 200, </>(x; v) — — In v + (f>o(x) + 
0{y) w 6.0, and ft 2 <T; v) = — (f)o(x) + ln2 ss —1.10. 


to the total drop /^(O) — ^2(1) = 0(—Ini?) 1 for small R. The only way to realize 
the large total drop p 2(0) — /^(l) is that p 2 (x; v) drops 0(— Ini?) over the subinterval 
(6,1). Indeed, this is the case as shown below. □ 


5.2.3 Internal dynamics over the interval ( b , 1) 


Proposition 5.9. For x E (6, 1), 


4>{x\v) = <j>o{x) + <t>i (x)v + 0(v 2 ), c 2 (x\v) = c 2 o{x) + c 21 (x)v + O(X), 


where 


, , x . C 20 {x) c 2 i{x) 2Jn c 20 (a;) 

0o(a;) = —In ———, q>\{x) =-— + —— In 


R 


C 2 o(x) J20 


R 


C 20 (x) =R+^{H{ 1) - H( x)), c 21 (x) = Jn t — (g(l) - #(*)); 


in particular, /i 2 (r; v) = p 2 o(x) + P 2 i(x)v + 0(z/ 2 ) where 


2c 2 i(x) _ 2Jn c 20 (x) 

C2o(^) ^20 i? 


M20 (#) = 2 In c 2 o (#) - In i?, /i 2 i (#) 
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Concentration of Anion in (b,1) Electropotential in (b,1) 



Electrochemical Potential of Anion in (b, 1) 



Fig. 8. Profiles of C2(x;is), (j)(x;u), and g2(x;n) over interval [6,1]. For v — 10 2 and 
R = 10 -5 small used for numerics, C2(x\v) = 0(\fR,v), <j)(x\v) = <j>o(x) + 0(y) ~ —5.0, 
and ^2(6; v) ~ —1.10 and /X2 (1) = In R « —11.5. 


The profiles of c 2 (x\u), 0(x; z^), and 1 x 2 (x; over (6,1) are shown in Figure 8. 

Corollary 5.10. Over the interval (6, 1 ), for small v and R, c 2 q(x) changes from 

c 2 o(fr) = 0(VR) to C 2 o(l) = R monotonically, and Ji 2 (x\v) changes from Jl 2 (6; z/) = 
0(1) to jd 2 ( 1 ) = lnR. Therefore, for x E (6,1), from (20), 

J 2 « 0(V~R) and jd 2 (x\v) drops O(-lnR) over (6,1). 

Remark 5.11. Note that, over this interval, the order of jd' 2 (x]v) varies in x from 

jd' 2 (b; v) = 0(1) to jd' 2 (l; u) = 0(1/V~R) but overall drops is 0(—In 17). This is different 
from what happened over the intervals (0, a) and (a, 6). □ 


5.3 Summary of mechanism for declining phenomenon 

In summary, for small v and R , with the technical assumption that v In R < \[R, we 
have, J 2 = 0(VR) over the whole interval (0,1) but with completely DIFFERENT 
scenarios over different subintervals (0,a), (a, 6) and (6,1). More precisely, 
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Fig. 9. Profile C 2 <T; v) over [0,1]. 


Concentration of Anion in (0,1) 


(i) over (0,a), one has C 2 (x;is) = 0(1) but fi' 2 (x',v) = O(VR) C 1 so that J 2 = 

0(y/R); (Note that the drop of ]±2 over the interval (0,a) is of order 0(VR), 
which has nearly no contribution to the drop of ]±2 over the whole interval (0,1).) 


(ii) over (a, 6), C 2 (x',u) = 0(1/u) but p' 2 (x',u) = 0(o\/R) so that J 2 = 0(VR ); 

(Note that the drop of JI 2 over the interval (a, b ) is of order 0(v\/~R), which is 
even smaller than that over the subinterval (0, a ) and, of course, has nearly no 
contribution to the drop of jd 2 over the whole interval (0,1).) 


(iii) over ( b , 1 ), different from what happened over each of the previous two subin¬ 
tervals, the orders of C 2 (x ; v) and jd 2 (x; v) are NOT uniform for x E ( b , 1 ) but 
the drop of O(lni^) in p2(x; v) is fully realized over this subinterval (b , 1) (see 
Remark 5.11). 

We comment that, for large Q 0 > 0, when R is small, both the transmembrane 
electrochemical potentials jd i(0) — jdi(l) = V + InL — In R for the cation and /^(O) — 
jd 2 ( 1 ) = — V -T In L — In R for the anion are large. The drop of pi(x‘,v) is mainly 
realized over the interval (a, b) but the drop of p2(x m ,u) is mainly realized over the 
interval ( 6 , 1 ). 

Here we provide the profiles of concentration (Fig. 9), electric field (Fig. 10) and 
electrochemical potential (Fig. 11) of the anion over the whole interval [0,1]. The 
functions C 2 (x ; v) and <f(x; v) over interval ( 0 , 1 ) are not continuous, because we make 
the plots of system (1) with 5 = 0. For the limiting system, there are fast layers at 
x = a and x = b where C 2 (x ; v) and <f(x; v) are discontinuous, but ft 2 (x; v) keeps the 
same value in fast layers. 






mu2 M phi(x) 
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Electropotential in (0,1) 



Fig. 10. Profile of <j)(x\ v) over [0,1]. 

Electrochemical Potential of Anion in (0,1) 



Fig. 11. Profile of /i 2 <A; v) over [0,1]. 


6 Concluding remarks 

In this work, we examine effects of large permanent charges on ionic flow through ion 
channels based on a quasi-one dimensional Poisson-Nernst-Planck model. We show 
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that one of the defining properties of transporters, obligatory exchange, can arise in 
an open channel with just one structure. When the permanent charge is large, the 
current carried by counter ions, majority charge carriers with the opposite sign from 
the permanent charge, can decline, even to zero, as the driving force (the gradient 
of electrochemical potential) increases. We also show that large permanent charges 
essentially inhibit the flux of co-ions, regardless of the magnitude of transmembrane 
electrochemical potential. 
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